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Abstract 

Skyrme theory is among the viable effective theories which emerge from low- 
energy limit of quantum chromodynamics. Many of its generalizations include 
also a dilaton. Here we find new self-gravitating solutions, both solitons and 
black holes, in a Generalized Skyrme Model (GSM) in which a dilaton is present. 
The investigation of the properties of the solutions is done numerically. We find 
that the introduction of the dilaton in the theory does not change the picture 
qualitatively, only quantitatively. The model considered here has one free pa- 
rameter more than the Einstein-Skyrme model which comes from the potential of 
the dilaton. We have applied also the turning point method to establish that one 
of the black-hole branches of solutions is unstable. The turning point method 
here is based on the first law of black-hole thermodynamics a detailed derivation 
of which is given in the Appendix of the paper. 
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1 Introduction 



Quantum chromodynamics (QCD) is a very successful theory in the regime of high 
energies where the relevant degrees of freedom are quarks and gluons. This is the so- 
called asymptotic freedom phase. In the low-energy regime, however, QCD is extremely 
hard to treat since in that regime it is nonperturbative. As it was shown by t'Hooft 
[T] QCD can be treated perturbatively in the approximation of large color number 
Nc -T- oo. Within that approximation QCD is equivalent to a local field theory of 
mesons and "glueballs" . Latter, Witten |2l |3] showed that baryons can be represented 
by topological solutions - solitons, in an effective theory of mesons. Such an effective 
theory had actually been proposed much earlier by Skyrme |H|5]. In his theory baryons 
emerge as topological solitons, the so-called skyrmions. Their conserved topological 
charge, the winding number, is identified as the baryon number. 

Skyrme theory has proved successful in the description of the static properties of 
baryons but it has also many deficiencies. Among them are the large masses and the 
missing intermediate-range nucleon-nucleon attraction which is vital for the formation 
of nuclei. Different generalizations of the theory have been suggested to remedy these 
deficiencies. One of the approaches is to add higher order derivative terms with the 
proper sign in the Lagrangian. Another approach is the inclusion of additional fields in 
the theory. The scale invariance and the trace anomaly of the underlying QCD theory 
in the low-energy limit leave place for the inclusion of scalar particles. For that reason 
many of the suggested generalizations of Skyrme theory include a dilaton. The dilaton 
has considerable contribution to the intermediate-range attractive forces [HI H El IS]- 
The Lagrangian of a Generalized Skyrme Model (GSM) including a dilaton has been 
derived from QCD in the low-energy regime in [lOl [111 112] . 

Skyrme theory coupled to gravity to our knowledge has been considered for the first 
time by Luckock and Moss [13]. They studied black-hole solutions. Self-gravitating 
skyrmions have been studied by Glendenning, Kodama, and Klinkhamer [Hj in the ap- 
proximation of large baryon numbers. They found that gravitating skyrmions cannot 
be applied as models of baryon stars since they are energetically unstable and disinte- 
grate into separate nonbounded solitons. Self-gravitating solitons with baryon number 
one have been investigated numerically by Droz, Heuseler, and Straumann p5]. They 
have studied not only regular solutions but also black holes. Their motivation was to 
find stable black-hole solutions which serve as counterexamples of the no-hair conjec- 
ture. The black holes that they obtained have nontrivial Skyrme hair. The energy 
density of the matter field (the Skyrme field) decays rapidly with the increase of the 
distance between the observer and the source and a distant observer cannot distin- 
guish them from the Schwarzschild black hole. In other words the Skyrme field does 
not contribute to new independent, asymptotic charges. The authors found that these 
solutions are stable against spherically sjTiimetric perturbations [161 lH]- Soon after 
that, Bizon and Chmaj [18] found another branch of solutions in the same model which 
turned out to be unstable. The stability of Einstein-Skyrme (ES) solitons and black 
holes has been studied also in [THl 1201 EI]- The Einstein-Skyrme black holes have been 
generalized to the case of asymptotically anti-de Sitter spacetime in [221 123]- Solitons 
and black holes in SU{3) ES theory can be found in [24j. The interior of the ES black 
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holes has been studied in [25]. 



All of the gravitating solutions discussed so far are spherically symmetric and are 
based on the so-called hedgehog ansatz. Monopole black-hole skyrmions based on a 
nontopological ansatz for the Skyrme field have been constructed in [26] . Axisymmetric 
soliton and black- hole solutions are given in [27] . The so-called harmonic map ansatz, 
which is a generalization of the rational map ansatz, has been applied in [28] to obtain 
axially symmetric SU{3) gravitating skyrmions and in [29] to obtain self-gravitating 
SU{2) skyrmions with discrete symmetries - platonic skyrmions. The same authors 
obtained also solutions for spinning skyrmions [30]. Employing the rational map ansatz 
Piette and Probert constructed self-gravitating configurations which, unlike the solu- 
tions obtained with the hedgehog ansatz [H], remain bounded even for very large 
baryon numbers [3ljl|. 




Self-gravitating objects in the generalizations of the Skyrme model which include a 
dilaton have not been so intensely studied. Kalbermann obtained an effective equation 
of state (EOS) for nuclear matter starting from a generalized Skyrme Lagrangian in 
which a dilaton and a co meson are included [32] [§. In this approach nuclear matter is 
described as skyrmion fluid which consists of free skyrmions with baryon number one 
immersed in a mean field background. With this EOS massive skyrmion stars were 
constructed in [33]-[2Z] Ifl and they were compared to neutron stars based on modern 
equations of state. Astrophysical applications of the solutions and present plausible 
skyrmion star candidates are also given there. It appears that skyrmion stars are quite 
promising since they can explain the high masses of some stars that have been observed 
and cannot be explained by some of the modern EOS. Massive stars in theories with 
stiff EOS, such as Skyrme's, have been discussed from an observational perspective in 
[38]. 

To our knowledge, self-gravitating soliton and black-hole solutions coupled to a 
dilaton have not been studied so far. The purpose of the present paper is to make a 
step in that direction. We generalize the paper of Bizon and Chmaj [18] by adding 
a dilaton in the model. In particular we study self-gravitating solutions in the GSM 
proposed in [12], [39] . 

The paper is organized as follows. The model is presented briefly in Sec. [2l The 
results of the numerical investigation of the solutions are presented in Sec. [3] in which 
subsection 13.11 is dedicated to the solitons, and subsection 13.21 - to the black-hole 
solutions. In Sec. ID the stability of the GSM black-hole solutions is studied with the 
application of Poincare's turning point method. Here the turning point method is 
based on the first law of black-hole thermodynamics a detailed derivation of which is 
given in the Appendix. The results are summarized in the Conclusion. 

-'^ Another possibility for the stabihzation of the muhisohton stars is the introduction of a chemical 
potential [40] , 

^The Lo meson has been considered also in [HI for the stabilization of the skyrmion without 
higher derivative terms. 

•^A wonderful brief review on Skyrme theory and its applications can be found in [35] 
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2 Basic equations 



We start with the following action: 



s= I d'^V^ai-^ + L,,]. (1) 



For the matter sector Lm we consider the effective GSM Lagrangian [12] in which the 
derivatives have been substituted with covariant derivatives 

1 TV f ^ 

Lm = -fl exp(-2a) TriV^f/V^fZ+l + exp(-2(T) g^^d^ad.a (2) 

+ ^Tr[(V,f/)f/+, {V^U)U^f + V{a). 

Here U is the SU(2) chiral field, a is the dilaton, is the covariant derivative with 
respect to the metric g^y, f-,, is the pion decay constant, e is the Skyrme constant, Cg 
is the gluon condensate, Nf is the number of flavors and e = 8Nf/{33 — "^Nf). 

The first two terms in ([2]) are the kinetic terms for the chiral and the dilaton fields. 
The third term is the one introduced by Skyrme for the stabilization of the soliton 
solutions and the last one is the potential of the dilaton field given by 
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exp(— 4cr) — 1 H — (1 — exp(— ecr)) 



(3) 



As mentioned above the dilaton is introduced in the theory to restore scale invariance. 
It couples to those terms of the Lagrangian density that break the scale invariance 
[U |l3] . In the classical Skyrme theory the first term breaks scale invariance while the 
stabilizing term is scale invariant. 

We will restrict our study to static, spherically symmetric, asymptotically fiat so- 
lutions. Hence, the following ansatz for the metric can be used: 

ds^ = A\r) ( 1 - de-(l- ' dr^ - r\de' + sin^ Ody^') (4) 



and the hedgehog ansatz for the chiral field 

U = exp[T-vF{r)], (5) 

where r are the Pauli matrices and f is a unit radial vector. It is also useful to introduce 
a new function $ for the scalar field instead of cr, defined by 

$ = exp(-(T). (6) 
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The reduced field equations following from the action ([T]) are 



where 



m' = -^^m, (7) 

A' = ^(^^F'^ + arN^'^'^A, (8) 

if^^n'-lfAF'^l^-lA^^O, (9) 

aiV(.VA.')'-^/A^-||-iA^ + ^M|K.o, (10) 



Cm = -f^F'' - 4 - aiV/$ '2 + V^($), (11) 



y($) = 167rG'V($) = -7iV 



$4_i + l(i_$^) 



(12) 



/ = 1-^^, M = ar2$2 + 26sin2F, = f 2 a + 6 ^in' i^- (13) 
r \ / 

We have introduced the following constants: 

a = MGfl, b = 8nG\, 7 = 2ttG^, N = (14) 



where a and are dimensionless. It can be easily seen that the system of field equations 
(17|)-( fT0|) does not admit in general solutions with trivial dilaton field a = ($ = 1). 

For the numerical treatment of the system it is convenient to work with the following 
dimensionless variables and parameters: 



X = ef^r, fi{x) = ef^m, D^s = 7rz~J7^- (^5) 



7A^ 
2^/2- 



The dimensionless field equations can be written in a form which depends only on the 
two coupling parameters a and D^s and on the parameter connected to the number 
of the flavors A^. For the number of flavors we fixed the value Nf = 2, so iV = 1. 

We will consider two classes of solution to the above field equations - solitons and 
black holes that can be obtained by imposing the appropriate boundary conditions. In 
the next section we will describe in detail the required boundary conditions and the 
obtained results. 



^The parameter a is two times bigger than the parameter a used in [18| and [15| and the parameter 
Doff is chosen to be the same as in p2] . |39j . 
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3 Numerical results 



3.1 Solitons 

The domain of integration for the sohton solutions is x G [0, oo) and the requirement 
that the functions should be regular at the origin leads to the following expansions of 
the functions at x ^ 0: 



aDefr($o'-l + -(l-'^o')l + 



4 
-I 

e 



(16) 



Fix) 



\x->-0 



Ao + \aAoF^' + 2F^') x' + 0{x'), 
nil + FqX + 0{x^), 



2D 



3N 



2N " 



x^ + 0(x^ 



(17) 
(18) 

(19) 



where n is an integer number, Aq and $o are the values of the corresponding functions 
at x = and Fq is the first derivative of F with respect to x at x = 0. The shooting 
parameters are Aq, $0 and Fq. 

The asymptotic flatness of the solutions requires [f| 



A{x)l 



F(x)L 



0. (20) 



The last boundary condition for the function F corresponds to the following require- 
ment for the Skyrme chiral field: U{x)\x-^oo = 1- This effectively compactifies to 
and the chiral field can be considered as a map from to SU(2), i.e., a map from 
to since SU(2) is topologically S'^. The third homotopy class of SU(2) is Z which 
means that every field configuration is characterized by an integer winding number B. 
This winding number is usually interpreted as the baryon number and is given by 



B 



where the topological current B^ is 



cfx^/^Bo, 



B„ 



247r2 



Tr[(f/+V''t/)(f/+V"t/)(f/+V^f/)] 



(21) 



(22) 



For the hedgehog ansatz and for the imposed boundary conditions B 



n. 



^The boundary conditions at infinity are also manifested by the asymptotic expansion of the 
functions at infinity (j37|) and 
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Figure 1: The metric functions fi{x) (left panel) and A{x) (right panel) for soliton 
solutions with different value of the parameter n where D^s = 0.00025, a = 0.15 and 
= 1. Solitons which belong only to the lower branch of solutions are presented and 
the upper branch of solutions have similar qualitative behavior {lower and upper branch 
are defined below). 



The field equations fl71)- f|T0l) together with the above boundary conditions are solved 
numerically using the shooting method. On Figs. [1] and [2] the metric functions fi{x) 
and A{x), the Skrymion field F{x) and the scalar field $(x) are presented for several 
soliton solutions obtained for Des = 0.00025, a = 0.15, = 1 and for different values 
of the parameter n. As it can be seen the profiles of the functions differ significantly 
for different values of n and a general observation is that the number of the "steps" of 
the function F{x) is equal to n — 1 

Similarly to the case when the dilaton is not present, nonuniqueness of the solutions 
exits. For fixed values of the input parameters two soliton solutions with different 
values of the shooting parameters exist which belong to the so-called upper and lower 
branches of solutionqj. These branches are presented on Fig. [3] where the shooting 
parameters and $o are given as a function of the coupling parameter a for several 
values of Dgff- The F^^a) and ^o{a) dependences for different values of n are given on 
Fig. a for Deff = 0.00025. 

The obtained sequences of solutions are qualitatively very similar to the case with- 
out dilaton field. The main difference is that the maximal value of the coupling param- 
eter Omax for which soliton solutions exist depends strongly on the coupling parameter 
Des and can be several times larger than in the case without scalar field. An interesting 
observation is that when Des approaches the values considered in [I2],[39] (on the plots 
Deff = 0.1), ttmax approaches the values in the Einstein-Skyrme model (ESM). 

^Such "steps" of the F{x) function exist also in the case when no scalar field is present. 

^When we refer to upper and lower branch of the soliton solutions, we use the convention from |18) 
which is based on the Fq phase diagram. This is exactly opposite to the $o phase diagram where the 
upper branch is below the lower branch. 
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Figure 2: The skyrmion field F{x) (left panel) and the dilaton (right panel) for 
the same solutions as on Figure [1] are shown. 




Figure 3: The shooting parameters and $o as functions of the coupling parameter 
a for different D^s are shown. 
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3.2 Black Holes 



In the case of black holes the domain of integration is x G [xh, oo), where xh is the 
dimensionless radius of the black-hole horizon. The presence of a horizon makes the 
black hole topologically trivial because the three-dimensional space is equivalent to 

minus a ball and the map U{x) is topologically trivial. That is why the function F 
does not have to be a multiple of vr at the black-hole horizon. 

The boundary conditions at infinity are again deduced from the asymptotic flatness 
requirement and they are the same as in the soliton case (120|) . The boundary condition 
for the dimensionless metric function fi at the horizon is the standard one 



The regularization conditions at the horizon are 



dF 

dx 



sin2F 



$2 



sin^ F 



X 



$2 



)sinj_F 



(23) 



(24) 



X 



1 - f sin^ F (^2$2 + 



sin^ F 

x"^ 



X- 



dx 



¥Xh 



Ax^{^' 
N 



£-l^ 



2 sin^ F 



N 



X 



(25) 



1 - f sin^ F ( 2$2 + 



sin'' F 



- x^aD^s - 1 



4(1 



Using the above boundary conditions we can conclude that in the case of the black 
holes we have again three shooting parameters, A^, F^, and which denote the 
values of the corresponding functions on the horizon. 
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Like in the ESM [18], in the GSM two branches of black-hole solutions exist which 
we will denote by upper and lower brancl^. These branches are presented on Figs. 
El and E] where the shooting parameters and are plotted as functions of the 
radius of the horizon for sequences of black-hole solution when we vary the coupling 
parameters Dcs and a. The metric functions fi{x) and A{x), the Skrymion field F{x), 
and the scalar field $(a;) are presented on Figs. [7] and [8] for several black-hole solutions 
of the lower branch with a = 0.15, = 0.05 and for different -Des- The qualitative 
behavior of these functions is the same for the upper branch. 

As it can be seen for fixed values of the coupling parameters, black-hole solutions 
exist up to a certain value of the radius of the horizon x^"^. For fixed values of 
Des, the maximal radius of the horizon x^"^ decreases with the increase of a and 
eventually reaches zero at some Cmax which means that black-hole solutions do not 
exist for a > flmax- The value of Omax for which x^""^ goes to zero decreases when we 
increase Defi- 

Like in the soliton case, the main difference between the black-hole solutions with 
and without scalar field is not qualitative but quantitative. The values of x'J}"'^ and 
Omax depend strongly on the parameter D^s and they can be several times larger in the 
GSM than in the ESM. Again if we consider D^s close to the one used in [12], [31] the 
values of x^""^ and Omax approach the corresponding values in the ESM. 

There are some interesting properties of the solutions in GSM which are also present 
in ESM. For fixed D^s, the Omax for black holes is the same as Omax for the n = 1 solitons 
within the numerical error. When the radius of the horizon goes to zero xh the 
masses of the black holes M and the values of the shooting parameters and 
approach the values of M, Fq and $o of the corresponding n = 1 solitons. The black- 
hole solutions presented so far correspond to the ground state (baryon number one) 
solitons. An infinite series of excitations of both the upper and the lower branches of 
solutions corresponding to the solitons with > 1 is also observed. As an example the 
first excitation is shown on Fig. M (the dashed line). As it is expected, the values of 
M, Fh and approach the corresponding values of the n = 2 solitons when xh — ^ 0. 
The results suggest that an infinite number of excited black-hole solutions exists and 
they are in one to one correspondence with the excited soliton solutions with n > 1. 

4 Stability and thermodynamics of the black holes 

A very convenient method to infer information for the presence of instabilities in the 
black-hole solutions is the Poincare's turning point method [44j. We will apply it here to 
prove that the lower branch of black-hole solutions is unstable. The method is based on 
the thermodynamics of the black holes. According to that method a change of stability 
is indicated by a turning point or a bifurcation point on the proper conjugate diagram. 
In the vicinity of a turning point the branch with a negative slope is unstable. Further 

^The notations upper and lower branch are chosen using the Fh{xh) dependences. This is obviously 
different to the soliton case and actually the upper branch for black holes will have similar properties 
(such as stability, finiteness/divergency of some of the functions as a — > ) to the lower branch for 
solitons and vice versa. 



10 




Figure 5: The shooting parameter Fh as a function of the radius of the black-hole 
horizon xh for several values of the coupling parameters. The results for Dgfr = 0.005 
and for various values of a are shown on the left panel and the results for a = 0.1 and 
for various valued of the parameter Des are shown on the right panel. 




Figure 6: The shooting parameter as a function of the radius of the black- hole 
horizon xh for the same black-hole solutions as on Fig. [5] is shown. 
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Figure 7: The metric functions (left panel) and A{x) (right panel) for black-hole 
solutions with different values of Des (a = 0.15 and = 1). Solutions belonging only 
to the lower branch of solutions are presented since the upper branch of solutions have 
similar qualitative behavior. 
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2-1 . 1 . 1 . 1 . 1 . 1 0.00 -F 1 I . 1 — . 1 . 1 . 1 

0.00 0.05 0.10 0.15 0.20 0.25 0.00 0.05 0.10 0.15 0.20 0.25 

Figure 9: The shooting parameters Fh and $h as functions of the radius of the black- 
hole horizon xh for D^.^ = 0.00025 and a = 0.1. Two sequences of black-hole solutions 
are shown where the thick line corresponds to the same black-hole solutions as the 
thick line on the right panels of Figs. E] and El and the dashed line shows the excited 
black-hole solutions corresponding to the n = 2 soliton solutions. 



details both on the formal proof and on the application of the method for the study of 
the stability of compact objects in gravity can be found in |13]-[Sn]- Some aspects of 
the thermodynamics of the Einstein- Skyrme black holes have been considered in [5Tj . 

The first law of thermodynamics (FLTD) for the Einstein-Skyrme black holes has 
been derived in [52]. A more general derivation has been given later in [53]. For the 
black holes studied here the FLTD has the form 

dM = T dS. (26) 

A detailed derivation of (!26|) is given in Appendix |Xl In a micro-canonical ensemble the 
equilibrium solutions are the extrema of the entropy dS = dM, so the conjugate 
variables are the mass M and the inverse temperature T^^. The conjugate diagram 
T~^(M) is given in Fig. [10] for the black-hole solutions presented on Figs. [5] and [6] 
(right panels). It can be seen that the lower branch and the upper branch merge at 
a turning point. The slope of the lower branch on the conjugate diagram is negative 
which, according to the turning point method, means that the branch is unstable. It 
is expected that the excited black-hole solutions are dynamically unstable. 

The linear stability of both the soliton and the black-hole solutions presented here 
is currently being investigated. 

5 Conclusion 

In the current paper we have studied new self-gravitating soliton and black-hole solu- 
tions in a GSM. The numerical investigation shows that the inclusion of the dilaton 
does not change the qualitative picture - it remains the same as in the original ESM. 
There are some quantitative differences, though. 
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0.16 0.18 0.20 0.22 0.24 



M 

Figure 10: The inverse temperature as a function of the black-hole mass M for 
sequences of black-hole solutions for different D^s is shown. 



For both soliton and black holes nonuniqueness of solutions is observed. If we vary 
the parameter a and keep the rest of the parameters fixed we obtain two branches of 
solutions - the upper and lower branch. They merge at some maximum value Omax 
and no solutions are observed beyond it. The value of flmax of the solitons coincides 
with that of the black holes. The same coincidence is observed also in ESM. In GSM, 
however, amax is highly dependent on the additional parameter D^q which is related 
to the dilaton. It increases with the decrease of D^^ and may be several times higher 
than in the ESM. 

For the black holes, if we fix a and vary the radius of the event horizon xh the two 
branches merge at some maximal radius x'^'^^ . Beyond x^"^ no black-hole solutions 
exist. Again the quantitative difference between ESM and GSM is that due to the 
presence of the dilaton x'J^'^'-^ may be considerably higher in the GSM when the same 
values of a are considered. 

In both ESM and GSM with the decrease of the radius of the horizon r// — )■ the 
masses of the black holes approach the masses of the solitons from the corresponding 
branch of solutions. 

With the increase of the baryon number an infinite series of solitons and corre- 
sponding to them black holes exist which, however, are expected to be dynamically 
instable. 

In the present paper we also generalize the FLTD of black-hole thermodynamics 
for the GSM. It has the same form as in the ESM. From the application of the turning 
method we inferred information for the instability of the lower branch of black holes 
which is in agreement with the situation in ESM. 
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A Derivation of the First Law of Thermodynamics 



For the derivation of the first law of thermodynamics we follow the scheme in [52] . 
Following standard definitions for the Hawking temperature we obtain 



T = —f'A 

An 



1 



r=TH 



Attth 



1 



1 



Aim) [l-7;VH + 7;rHVH 



(27) 



which together with the entropy of the horizon S = irrfj gives 



T 



dS 
drn 



- A{rH) 1 - 7T + 7T 



H 



The mass of the black holes is given by the formula 

M = lim -r (1 - /). 

r— >oo 2 

In order to calculate the variation of the mass we will use the expression |^ 

rfA = j^ dr'A{r') (^1 - + ^r'^ , 



{2t 



(29) 



(30) 



which can be obtained with the help of eqs. and ([H]). Then the variation of M is 

(31) 



dM _ ^ as 

drn Oth 2 



I roo Q 

or I 



If the integral in this expression vanishes then 

dM ^ dS 
T- 



A{r') ll--v + -r'^V 



dr 



H 



dr 



H 



and the FLTD has the usual form 



dM = T dS. 



(32) 



The integrable expression can be brought to the form of full derivative with respect to 
the integration variable which allows the integral to be estimated. For the derivative 
in the integral the Leibnitz rule gives 



d 
drn 



A 1 --v + -r'V 



dA 



dr 



H 



1 V + - r^V \ + A — t:^ + - r — 



1 dv 



dV\ 9$ 



2 (9$ 



9$ / drH 



2 dF J dm 



(33) 



^Whcn r ^ oo, r J A = M + 0{r-'^). 
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We will treat each of the terms on the right-hand side of the expression above separately. 
From eq. ([H]) we obtain 



dA 



dr 



A 



H 



-F'- h aA^r $' 

r or 



H 



+ 



dr 



dr 



H 



dr 



H 



2r' dF 



uF' 



r 



+ 



+ 



^ 9$ 
dr - — 
orn 



J_F"^-aN{r'<^')' 
2r' 9$ ^ ^ 



For the expressions in the brackets of the second term and the third term on the 
right-hand side from the field equations for F and $ - ([9]) and ffTOj) . respectively, we 
have that 

-{fAuF')' + lfAF'^ 



1 dv 
2^dF 



and 



A - 



1 dv 



+ - 



dV 



2 9$ 2 9$ 
Collecting all terms for the integrable we obtain 



dF 

-aNir^ fA^')' + -fAF'^^. 



d 
drn 



A 1 



1 



1 



r^V 



(66)' 



where 
and 

6 = 



^^ = rfA = j^ dr'A(r') + 



(34) 
(35) 



dr 



, dF 
dru 



2r' dF 



uF 



+ 



+ 



dr 



, (9$ 

dr 



H 



J_F'^^-aN{r'<i>y 
2r' a$ ^ ' 



From we obtain for the integral 
d 



l^dr'^ A{r')(^l-^v + ^r"V^ = dr' {F,F^)' = F,F^ 



oo 



(36) 



^2 is regular on the event horizon so ^1^2 



0. The situation at infinity is more 



subtle since — )■ 00 and ^2|oo ~^ 0- The value of the integral at infinity can be 
estimated with the help of the following asymptotic expansion of the functions: 



m(r)| 



F(r) 



= M — 



+ 



1 — 



laFo 



2 r^ 



+ OKh(37) 



$(rl 



1 + -3 h O 

Njle - 4)r6 



(38) 



where F2 is a constant. With (1371) and fl38|) it can be found that 66 loo ~^ ^ which 
completes the proof of the FLTD fl32|) . 
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